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Abstract 

Polarized photon structure can be studied in the future polarized 
e + e~ colliding-beam experiments. We investigate the spin- dependent 
structure function of the virtual photon gj(x, Q 2 , P 2 ), in perturbative 
QCD for A 2 « P 2 « Q 2 , where — Q 2 (-P 2 ) is the mass squared of 
the probe (target) photon. The analysis is performed to next-to-leading 
order in QCD. We particularly emphasize the renormalization scheme in- 
dependence of the result. The non-leading corrections significantly modify 
the leading log result, in particular, at large x as well as at small x. We 
also discuss the non-vanishing first moment sum rule of gj, where 0(a s ) 
corrections are computed. 
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1 Introduction 



In recent years there has been growing interest in the study of a polarized photon 
structure function. The information on the spin structure of the photon would be 
provided by the resolved photon process in polarized electron and proton collision 
in the polarized version of HERA [l], || . More directly, the spin-dependent structure 
function of photon gj can be measured by the polarized e + e~ collision in the future 
linear colliders (Fig.l). 

From the theoretical viewpoint, the first moment of a photon structure function 
gl has recently attracted attention in the literatures || |], ^ ||, 0] in connection with 
its relevance for the axial anomaly, which has also played an important role in the 
QCD analysis of the spin structure of the nucleon. Our aim here is to carry out 
the QCD computation of photon's polarized structure function at the same level 
of unpolarized case. Here we note that the two-loop splitting functions of DGLAP 
equation or equivalently the two-loop anomalous dimensions have recently been 
calculated || |§, and we can perform the next-to- leading order QCD analysis for 
the polarized photon structure function. Actually there has already been an analysis 
of spin-dependent structure function gl for the real photon target by Stratmann and 
Vogelsang |10[ . 



In this paper we shall investigate the polarized virtual photon structure function 
gJ(x,Q 2 ,P 2 ) to the next-leading order (NLO) in QCD, in the kinematical region: 

where — Q 2 (-P 2 ) is the mass squared of the probe (target) photon, and A is the 
QCD scale parameter. We can base our arguments either on DGLAP-type Q 2 evolu- 
tion equation for the parton distributions or on the framework of operator product 
expansion (OPE) and the renormalization group (RG) method. The unpolarized 
virtual photon structure functions F%{x, Q 2 , P 2 ) and F2(x,Q 2 , P 2 ) were studied in 
the leading order (LO) |TT| and in the NLO ||T2| , |13f . And the parton contents of 



virtual photon were studied in ref . p4| , |T5f and the target mass effect of unpolarized 



and polarized virtual photon structure in LO was discussed in ref. |16| . 

The advantage to study the virtual photon target is that we can calculate the 
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whole structure function entirely up to NLO by the perturbative method. On the 



other hand, for the real photon target [p~T| , |T8|| we can calculate the perturbative 
pieces, but not the non-perturbative contributions which may be estimated, for 



example, by vector-dominance model fl9|| . The perturbative pieces for the real 
photon target can be reproduced from the result for the virtual photon case. 

In the next section we discuss the polarized photon structure functions. Next we 
present the two theoretical frameworks based on OPE (section 3) and on DGLAP 
parton model approach (section 4). In section 5, the sum rule for the first moment 
of gj will be evaluated up to the order of a s . The numerical analysis of gj will be 
given in section 6. The final section is devoted to the conclusion and discussion. 



2 Polarized photon structure functions 

Let us consider the forward virtual photon scattering amplitude (Fig. 2), 

T^rrfaq) = i J rf 4 xdV 4 ze^e^^(0|T(J M (a;)J ; ,(0)J p ( 2 /)J r (^))|0). (2.1) 

Its absorptive part is related to the structure tensor W^^lp, q) for the photon with 
mass squared p 2 = —P 2 probed by the photon with q 2 = —Q 2 : 

W^ up r{p,q) = -ImT^prtp, q) . (2.2) 

7T 

The anti-symmetric part, W^, pT , which is antisymmetric under the interchange of \i 
and z/, can be decomposed as 

+ t^x a q X {p ■ q epr^pp - t pTaP p p p° r q a ) 1 ' gl , (2.3) 

which gives two spin-dependent structure functions, gj(x, Q 2 , P 2 ) and g\ (x, Q 2 , P 2 ). 
For a real photon, g^ is identically zero, and there exists only one spin struture 
function, gJ(x,Q 2 ). On the other hand, for the off-shell or virtual photon (P 2 ^ 0) 
target, we have two spin- dependent structure functions gj and gl- More detailed 
argument on the structure functions is given in the Appendix D. The gl is related 
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to the structure function W2, which was discussed some years ago in p0| , pTf , such 



that gi(x,Q 2 ) = 2W^(x,Q 2 ). Here we note that the LO QCD correction was first 



studied by one of the authors in [^2j and later in [^3], [J. 

First we note that the same framework used in the analysis of the nucleon spin 
structure functions can be applied in our case. We can either base our argument on 
the OPE supplementd by the RG method, or on the DGLAP type parton evolution 
equations. It should be noted the next-to-leading order analysis is now possible since 
the two-loop anomalous dimensions of twist- 2 operators R % n relevant for the spin- 
dependent structure function (or equivalently two-loop parton splitting functions) 
were calculated independently by two groups, by Mertig-van Neerven || and by 
Vogelsang ||. 

3 Theoretical framework based on OPE 



In our previous paper [2^], we based our argument on the QCD improved parton 
model approach. Here we start with theoretical framework based on the OPE and 
RG method. Applying OPE for the product of two electromagnetic currents we get 



• / d 4 xe^T( J,(x) J U (0)) A = -ie» v x„q X £ ( Jj)<fa ■ ■ ■ £ ^XJ 

n=odd \*v J i 

-i(e w xaquq p -e upX aq^q p -q 2 e^\a) £ (oaW ' ' ' q ^ S C 2,n i ^ W ^""M 3 - 1 ) 

n=odd \ / i 

For polarized deep inelastic scattering, the twist- 2 and twist-3 operators: R™, R% 
contribute to the structure functions in the scaling limit. For g\ only twist- 2 oper- 
ators are relevant. Now we can write down the moment sum rule for g\: 

[ 1 dxx n - 1 gUx,Q\P 2 )= £ Cj(g a / A i 2 ^0* 2 ) l a)( 7 (p)|iJiO' 2 )l7(p)>, (3-2) 



i=ip,G,NS,-y 



where R l n and C % n are the twist- 2 operators and their coefficient functions (Hereafter 
we suppress the index 1 for twist-2 operators), with fi being the renormalization point 
and a = e 2 /47r, the QED coupling constant, if}, G, NS and 7 stand for singlet quark, 
gluon, non-singlet quark and photon, respectively. The relevant twist-2 operators 
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K (i = *P(S), G, NS, 7) are given by [[22|: 

^W-Mn-l = jn-lfyiajy* . . . ^-i}^ ^ _ trace termS; (3,3) 

j^i-t^-i = 1 ^-i e {° a/37 G^ 1 J D^--- J D^- l} G' /37 -trace terms (3.4) 
R w~nn-i = in-i^Wpm . . . d^-i^q^ _ ^ _ trace terms ( 3 5 ) 

i^w-^-i = i^e^^F^d^ 2 ■■■d^- l} F^ -trace terms (3.6) 

where { } means complete symmetrization over the Lorentz indices 07/1 ■ • • 
.D^ 1 denotes covariant derivative, 1 is an nj x rij unit matrix and Q 2 h is the square 
of the nf x nf quark-charge matrix, with rif being the number of flavors. Here we 
note that the essential feature in the above equation is the appearance of photon 
operators in addition to the hadronic operators. 

For — p 2 = P 2 ^> A 2 , we can calculate the photon matrix elements of the hadronic 
operators perturbatively. Choosing /1 2 to be close to P 2 , we get, to the lowest order, 



a ( 1 P 2 \ 

(j(p)\RMh(p)) = — (--i^ln— + Alj , i = i>(S),G,NS (3.7) 

where = (K^) 1 are one-loop anomalous dimension matrix elements between 
the photon and hadronic operators. On the other hand, in the leading order of the 
QED coupling constant, a, we have for the photon operator K^: 

(7(p)K(/i)| 7 (p)) = l (3.8) 

It should be noted that the finite term A\ depends on the renormalization scheme 
for the operators R l n . Putting /x 2 = — p 2 = P 2 , we have 

{i(jp)\rM\i{p)%*=p* = > (3-9) 



and the n-th moment with this choice fi 2 = P 2 in (|3.2| ) becomes 

/ dxx n - x g{{x,Q 2 ,P 2 ) (3.10) 
Jo 



E (l(p)\K(^ = P 2 )\l(p))[Tex P 



, in(g) 

dg —. . 

s(Q 2 ) P{g) 
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The evolution factor in the last equation is found to be [|H]j : 

Texp 



*( p2 ) in{g) 

dg —. . 
<?(Q 2 ) /%) 










1 — 1 



(3.11) 



where 



M n (Q 2 /P 2 ,s(P 2 )) 
X n (Q 2 /P 2 ,2(P 2 ),a) 



Texp 



stp \m 

S(Q 2 ) P{9) 



3(? 2 ) K n {g,a) m 
dg — —-t — 1 exp 

s(Q 2 ) 13(g) 



2 ) 



, in{g') 
1 M 



(3.12) 



with 7„ and K n the hadronic anomalous dimension matrix and the off-diagonal 
element representing the mixing between the photon and hadron operators (see 
Appendix A). Thus we get 

f'dxx^gUx^^P 2 ) = ^A n - M n (Q 2 /P 2 ,g(P 2 ))C n (l,g(Q 2 )) 

JO 47T 

+ X n (Q 2 /P 2 ,g(P 2 ),a) • C n (l,g(Q 2 )) + CI (3.13) 



with 

A n = (A n , A^, A n ) . 
The coefficient functions are given by (see Appendix C), 



(3.14) 



C n (l,g) = CC(l,g) 



X 9 Dn 



a: 



CZ(l,g,a) = —8 7 B, 



(3.15) 



with 5^ =< e 2 >= Yh=i e l/ n f,^NS = 1,S 7 = 3n f < e 4 >= 3YhIi A- 
We then derive the following formula for the moments: 



f 1 dxx^gJix^Q 2 ^ 2 ) 
Jo 

^„ a 1 l + \?/2(3 a s (Q* 



/o 
a 1 



4tt 2/3, 



«=+,-,A r S 



«s(g 2 
MP 2 



i=+,-,A r 5 [ \ a; s(-P 2 ) / 



i=+,-,NS 



M p2 ). 



.(3.16) 



Here a s (Q 2 ) = g 2 (Q 2 )/47r is the QCD running coupling constant. In ( |3.16| ), we have 
defined 

f? = K° n P?C n (l,0) (i = +,-,NS). (3.17) 

where P"'s are projection operators given in the Appendix A. The coefficients A™, 
B™ and C n are computed from the NLO perturbative calculation, and are given by 

An _ K \- Pfln^Pj 1 r ( , q\ 1 K Q Pi pnr 1 ~ A"/2/3p 

+<P*C n (l, 0)-^— - 2/3 A n Pf C n (l, 0) (3.18) 

pn-(l) pn -i 

m = kIY\ n ? 3 \ c n (i,o)- 



; 2/3 + A«-A« " v ' y l + Af/2/? 



+ ^" 1 i:IJ T ™ " ^"iS (319) 

= 2/3 (5 7 ^ + A n -C n (l,0)) (3.20) 

where A™ (z = +, — , ATS 1 ) are eigenvalues of the 1-loop anomalous dimension matrix 
7^°) and are given in the Appendix A. p and the (3\ are the one- and two-loop (3 
functions, and /3 = 11 — 2nj/3 and Pi = 102 — 38rij/3. 

All the quantities necessary to evaluate P™, A™ , B™, and C n are now known and 
will be presented in Appendix. Two-loop results || || have been calculated in the 
MS scheme [p5[| . Actually the expressions of Eqs.( [3.16| ) and (|3.17|) -( |3.20|) are the 



same in form as the ones obtained before by one of the authors and Walsh for the 
case of the virtual photon structure function F% |12] . The explicit expressions of the 



one-loop and two-loop anomalous dimensions || |9] as well as one-loop coefficient 
functions |26], |27|, P8"|, P9[ pO|, |8], |9|| are given in the Appendix B and C. 



Eq.( |3.16D is our main result of the present paper. The first term is the LO result, 
and the remaining terms are the NLO QCD corrections. 

Now let us examine the renormalization scheme independence of the coefficients; 
Af, B™ and C n . As in the unpolarized case, £>™ can be written in terms of a 
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scheme-independent combination of 2-loop anomalous dimensions and 1-loop co- 
efficient functions in the hadronic sector. Using the scheme-independent coefficients 
R 2 n (|31|, |32], we can write 



£r = A"4 n {i = +,-NS) (3.21) 
where the explicit form for R l 2n is given in Eqs.(9)-(12) of ref.^TJ (see also [|l2j) and 

L ' = ^TT^m (3 ' 22) 

which is the coefficient of the leading-log term. The scheme independence of £>™ 
follows from these two equations. 

Regarding C n , we first consider the photon matrix elements of the renormal- 
ized quark and gluon operators. The finite matrix elements A n and the tree-level 
coefficient functions C n (l,0) are given by 

A n = 6(< e 2 >, 0, < e 4 > - < e 2 > 2 )i^ G (3.23) 



C n (l,0) 



/<e 2 > 

| (3.24) 

V i 



Hence we have 
Noting that 

we find C n equal to be 



A„-C„(l,0) = 6<e 4 > Al G (3-25) 
= — B G , <L = 3n f < e 4 > (3.26) 

C" = 12ft < e 4 > (Sg + iJo) . (3.27) 



Since the combination B G + yi^ G is scheme independent p5 |, so is C n . In fact, in 



the MS scheme ||, the gluon matrix elements of quark operators (i = tp, NS) read 
(p,s|fl<"" (^)\p,s) = A(p 2 ,fi 2 ,9)[{s'pf"-P''"- 1 }- traces] (i = i>,NS) 
A * = ^(^"'"I?^-) < 3 ' 28) 



where the finite matrix element A^ G is given in the parton language as 
where 

a sl 9 ( x ) = T /[( 4 - 82) {in x + ln(l - x)) 
with Tf = rif/2. Thus we get 



Ate = 2n, 



n — 1 



SAn) + 



1 



_n(n + 1) " v ' ' (n + l) 2 n 2 n. 
Therefore, from ( |3.27|) we finally arrive at 



C n = 24(3 f < e 4 > 



2 4 
+ 



n n + 1 n 2 (n + 1) 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



which is consistent with the Box diagram calculation. 

On the other hand, in the RG scheme adopted by Kodaira [28] which is the 
momentum subtraction scheme, we have for n > 3, 



with the coefficient function given by 



Bl 



2n 



f 



2 

^7 + 



n n + 1 n 2 ' (n + l) 2 
For n = 1, due to the Adler-Bell-Jackiw anomaly 



.4 



n=l G 



-2n 



f 



and 



B n G =1 = , 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



by definition, since we have no gauge-invariant gluon operator for n = 1. Combining 
this with the result for A^ G , we arrive at the same result for C n . 

The scheme- independence of the remaining coefficients A™ follows from the above 



arguments on and C n and the physically measurable moments given in Eq.( 3.16 ). 
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4 QCD improved parton model approach 



We now turn to the analysis based on the QCD improved parton model |[4]] using 
the DGLAP parton evolution equations. 

Let q±(x, Q 2 , P 2 ), G±(x,Q 2 ,P 2 ), T].(x,Q 2 ,P 2 ) be quark with i-flavor, gluon, 
and photon distribution functions with ± helicities of the longitudinally polarized 



virtual photon with mass — P 2 [24] . Then the spin-dependent parton distributions 
are defined as 

Aq l = qi + ql-ql-ql (4.1) 

ag 7 = g\ - ai , Ar 7 = r 7 - ri . (4.2) 

In the leading order of the electromagnetic coupling constant, a = e 2 /4n, AT 1 does 
not evolute with Q 2 and is set to be Ar 7 (a;, Q 2 , P 2 ) = 5(1 — x). The quark and 
gluon distributions Aq % and AG 7 satisfy the following evolution equations: 

+ AP qG (-,Q 2 ) AG^y,Q 2 ,P 2 )} + AP qij (x,Q 2 ,P 2 ), (4.3) 



+ AP GG (-,Q 2 ) AG\y,Q 2 ,P 2 )} + AP G ,(x,Q 2 ,P 2 ),(4.4) 

y J 

where APab is a polarized splitting function of 5-parton to A-parton, defined as 
AP AB = Pa+b+ — Pa-B+ (= Pa-B- — Pa + b~, due to parity conservation in QCD 
and QED). 

For later convenience we use, instead of Aq l , the flavor singlet and non-singlet 
combinations defined as follows: 

Ag 7 = £ Aq* (4.5) 

i 

Aq^ = Atf-^ (4.6) 
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so that J2i.Aq l NS = and n/ is the number of relevant active quark flavors. The 
quark-quark splitting function AP q i g j is made up of two pieces, the one representing 
the case that j-quark splits into i-quark without through gluon, and the other one 
through gluon, and may be expressed as 

AP qiqj = S tj AP„ + -j- AP q s q (4.7) 

where the second term is representing the splitting through gluon, and AP qq and 
AP qq are both independent of quark flavor, % and j. It is noted that by construction 
AP qq is relevant for the evolution of flavor-siglet Aq$ and first appears in the order 
of a 2 . 

In the QCD improved parton model, which is based on the factorization theorem, 
the polarized virtual photon structure function gj(x, Q 2 , P 2 ) is expressed as 

gJ(x,Q 2 ,P 2 ) = f^- \T Aq l {y,Q 2 ,P 2 ) C l {-,Q 2 ) 

J* y l i y 

+AG\y, Q 2 , P 2 ) C G (-, Q 2 )} + C'(x, Q 2 ), (4.8) 

y J 

where C l , C G , and C 7 are the coefficient functions of 2-quark, gluon, and photon, 
respectively, and are independent of target photon mass P 2 . Up to one- loop level 
they are given by, 

C\z,Q 2 ) = e 2 {5{l-z) + ^±B q {z)} (4.9) 

C G (z,Q 2 ) = <e 2 >{0 + ^±B G (z)} (4.10) 

C 7 (z,Q 2 ) = ^n f <e 4 > BJz) (4.11) 
Air 

where < e 2 >= J2i e f/ n f an d < e 4 >= J2i e i/ n f- It is noted that B q (z) in Eq.( |4.9|) 
is independent of the quark flavor i. Since Ya Aq l C % is rewritten as 

i it 

= AqUy, Q\ P 2 ) <e 2 > [5(1 - X -) + ^P-B q (^)} 

+ E^ 2 A&sfa&^Ul --) + *@±B q {*j\, (4.12) 

I y 47r y J 
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we obtain 

gj(x, Q 2 , P 2 ) = f 1 ^ Uql{y, Q 2 , P 2 ) C S { X -, Q 2 ) + AG J (y, Q 2 , P 2 ) C G (*, Q 2 ) 
Jx y I y y 

+Aql s (y,Q 2 ,P 2 ) C NS (-,Q 2 )) + C\x,Q 2 ) (4.13) 

y J 

where we have defined 

C s (z,Q 2 ) = <e 2 >{<5(l-^) + 

a,(Q 2 



C JV6 (z,Q 2 ) EE 5(1 -z) + 



47T 



Ag 



NS 



2 A J 





(4.14) 




(4.15) 
(4.16) 



and B s {z) = B NS (z) = B q (z). From Eqs.(PlHPT) and (p^) , the evolution 
equations for Aq^, AG 1 , and AqJj S are now given by 



dAql{x,Q 2 ,P 2 " 
dlnQ 2 



x y 



, x 



AP qq (-, Q 2 ) + AP q s (-, Q 2 ) I AqUy, Q 2 , P 



l2 D 2\ 



+ n f AP qG (-,Q 2 ) AG^(y,Q 2 ,P 2 ) ] j + Y, Ap ^Q 2 



(4.17) 



rfAG^(x,Q 2 ,P 2 ^ 
cflnQ 2 



rfA^(^Q 2 ,P 2 ^ 
cilnQ 2 



^{ap^^a^Q 2 ,^ 
* y l y 



,x 



+ AP GG (- Q 2 ) AG\y, Q 2 , P 2 ) \ + AP G7 (x, Q 2 ), (4.18) 



^AP^-.g 2 ) Ag; 5 (y,Q 2 ,P 2 ) 

* y y 

+ J2e 2 {AP q ^(x,Q 2 ) - -Y,^P q ^,Q 2 )} ■ (4.19) 

i n f j 



Introducing a row vector Aq 1 = (Agjy, AG 1 , AqJf S ), the above evolution equations, 
Eqs.( [4.17D -( |4.19| ) are expressed in a compact matrix form. 

dAq 2z% ,p2) = Afc(x ' q2) + r ^ p2 ) Ap (-' q 2 ) ( 4 - 2 °) 



where the elements of a row vector Ak = (AKs, AK G , AKns) are 

AK S = E A ^7 ? AK G = AP Gl 

i 

AK NS = £ e 2 {AP^--]TAP g U ( 4 - 21 ) 

i U f 3 

Since AP ? ; 7 is proportional to e 2 , it is easily seen that AKs an d AK NS have factors 
rif < e 2 > and rif(< e 4 > — < e 2 > 2 ), respectively. The 3x3 matrix AP(z, Q 2 ) is 
written as 

(AP q s q (z,Q 2 ) AP Gq (z,Q 2 ) \ 

AP(z,Q 2 )= \AP qG (z,Q 2 ) AP GG (z,Q 2 ) (4.22) 

V AP^^Q 2 )) 

where 

AP q S q EE AP^ + AP*, AP qG EE f AP qG 

AP Gq ee AP Gq , AP GG = AP GG , AP q N q s = AP qq (4.23) 

Once we get the information on the coefficient functions in Eqs.( |4.9| )-( |4.11| ) 
and parton splitting functions in Eqs.( |4.21] )-( |4.23| ), we can predict the behavior of 
gJ(x,Q 2 ,P 2 ) in QCD. The NLO analysis is now possible since the spin-dependent 
one-loop coefficient functions and two-loop parton splitting functions are avail- 
able |8|, ||. There are two methods to obtain gj(x, Q 2 , P 2 ) in NLO. In one method, 
we use the parton splitting functions up to two-loop level and we solve numerically 
Aq 1 (x , Q 2 , P 2 ) in Eq. Q4.20D by iteration, starting from the initial quark and gluon 
distributions of the virtual photon at Q 2 = P 2 . The interesting point of studying 
the virtual photon with mass — P 2 is that when P 2 ^> A 2 , the initial parton distri- 
butions of the photon are completely known up to the one-loop level in QCD. Then 
inserting the solved Aq 7 (x,Q 2 ,P 2 ) into Eq. (|4.13|) , and together with the known 
one-loop coefficient functions we can predict gi(x,Q 2 ,P 2 ) in NLO. 

The other method, which is more common than the former, is by making use of 
the inverse Mellin transformation. From now on we follow the latter method. First 
we take the Mellin moments of Eq. ([4.13| ), 

f 1 dxx n ~ l gl(x, Q 2 , P 2 ) = Aq>, Q 2 , P 2 ) ■ C(n, Q 2 ) + C\n, Q 2 ) , (4.24) 
Jo 
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where we have defined the moments of an arbitrary function f(x) as 



f(n) 



dxx n -'f(x). 



(4.25) 



Comparing Eq. ( |4.24| ) with Eqs.(3.10) and (3.15), we can easily see the correspon- 
dence between the quantities in the QCD improved parton model and those in the 
framework of OPE as follows: 



Aq\n,Q 2 ,P 2 ' 



= E <7(P)|W = ^ 2 )|7(P)> Texp 

j=S,G,7VS, 7 V 



S{P2) d ^ 

s(Q 2 ) 9 P{g) 



3 1 



C(n,Q 2 ) = C n (lrg) 
C^(n,Q 2 )=C^l,g,a) 



{i = S,G,NS) (4.26) 
(4.27) 
(4.28) 



Henceforce we omit the obvious n-dependence for simplicity. We expand the 
splitting functions Ak{Q 2 ) and AP(Q 2 ) in powers of the QCD and QED coupling 
constants as 

aa s (Q 2 ) 



Ak(Q 2 ) 
AP(Q 2 ) 



a 



2n + (2tt)2 



-Afc« 



a s (Q 2 
2tt 



a s (Q 2 
2tt 



AP(D + . . . 



(4.29) 
(4.30) 



and introduce t instead of Q 2 as the evolution variable 



(4.31) 



Pa « S (Q 2 ) 

Then, taking the Mellin moments of the both sides in Eq.( f4.20|) , we find that 
Aq 7 (t)(= Aqr 7 (n, Q 2 , P 2 )) satisfies the following inhomogenious differential equa- 
tion isa M: 



dt 2?r a, L 



2Po 



0(a s 



+ 



A<f (t){AP(°) + ^[AP« _ A A P(°) 
{ 2-k 1 2f3 



0(a 2 )\ (4.32) 
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where we have used the fact the QCD effective coupling constant a s (Q 2 ) satisfies 
da s (Q 2 ) „ a s (Q 2 ) 2 q a s (Q 2 ) 3 . 

"dlnQ^ " _ A_ (4^)5" + " ' (433) 

with fa = 11 - 2n//3 and ft = 102 - 3871/73. Note that the P 2 dependence of Aq 1 
solely comes from the initial condition (or boundary condition) as we will see below. 
We look for the solution Aqr 7 (t) in the following form: 

A<f(i) = Aq^°\t) + A(f (1) (t) (4.34) 

where the first (second) term represents the solution in the LO (NLO). First we 
discuss about the initial conditions of Aqr 7 . 

In Section 3, we have observed that for — p 2 = P 2 ^> A 2 the photon matrix 
elements of the hadronic operators R l n (i = 4>(S), G, NS) can be calculated pertur- 
batively. Choosing the square of the renormalization point fi 2 to be close to P 2 , we 
obtained, to the lowest order 

<7(P) I KM I 7(P)> = ^ ir\«^ + O , i = i>(S), G, NS . (4.35) 

The i^°'*-terms and A^-terms represent the operator mixing between the hadronic 
operators and photon operators in the LO and NLO, respectively. The operator mix- 
ing implies that there exists quark distribution in the photon. When we renormalize 
the photon matrix elements of the hadronic operators at fi 2 = P 2 , we obtain 

(7(P)IWI7(P))U=P-^ (4-36) 

which shows that, at fi 2 = P 2 , quark distribution exists in the photon, not in the 
LO but in the NLO. Thus we have 

A^°)(0) = 0, A<f«(0) = -^4 n . (4.37) 

Explicit expressions of A n in the MS scheme are given in Section 3. 

With these initial conditions, we obtain for the solution Aq 7 (t) of Eq. ( j4.32| ), 

2Ap(°± \ 

' (4.38) 



Aq^(t) = -i^-a 
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A<f (1) (t) 



2 a j/Ve^-^ [AP« 

„ . „rni 

a s (t) 



A 

2(3 C 



+b 1 



AP (o) 



a s (t) 



-AP<°)r AP(°)t 



,(0) 



« s (0) 



where 



-Afc<°>- 



1 - 



2AP<°) 



A) 



b = ^ 



a: 
2^ 



2/?o 



+ 2 a 



A p(D - A A p(°) 

2/3 



J AP(°) 



(4.39) 



(4.40) 



(4.41) 



It is noted that the parton distributions Ag 7 (t) c?o depend on the initial conditions 
Aq 7 (0) = (a/47i)A n , but we have seen in Section 3 that the structure function 
gl(x, Q 2 , P 2 ) itself is independent of A<f(0) in NLO in QCD. 

The moments of the splitting functions are related to the anomalous dimensions 
of operators as follows: 



AP (o) 
Afc<°> 



--7° 
K° 

- £V ni 



A p(D 

Afc« 



o / n 



(4.42) 
(4.43) 



The evaluation of Aq 7{0) (t) and Aqr 7(1) (t) in Eqs. flPsD and (^39|) can be easily done 
by introducing the projection operators P/ 1 such as 



AP (o) 



pn pn 

i 3 



7 E A r ^ r 

^ i=+ — NS 



+,-,NS 



i^j, 



En 1 



(4.44) 
(4.45) 



where A" are the eigenvalues of the matrix 7°. Then the solution Aq" ( (t) of Eq. fl4.32p 
in the NLO is written as 



a 



8ttA> 



-A k° y p™ - 



1 H — 



a s (t) 



1 1+7 



« s (0) 



+ Ik 1 y p. n — - — + ^-k° y p n (1 — - — 
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" n j-2[3 + A™ - A™ A?/2/? ^° A " ^ ^ 



x 1- 



<Xs(t) 

a s (0) 



A. 
2^0 



f y ^ ^ g 1 A y0 v pn K/Wo 
+ I n tT 2/? Q + A?-A? i + ^L f3 n V 1 + ^ 



2&j ' " * ' 2ft 

... ... . "2^ 

x 1 - 



-or a (0). 

+ 2/? A n (4.46) 

Now inserting the above solution of Aqr 7 (t) into Eq. ([4.24|) and together with 
the information on the coefficient functions in Eq. (3.15), we reproduce the same 
formula for the moments of g2(x,Q 2 ,P 2 ) given in Eqs.(3.16)-(3.20), as for the case 
of OPE approach in the NLO. 



5 Sum rule for the first moment of gj(x : Q 2 , P 

The polarized structure function gj of the real photon satisfies a remarkable sum 
rule |, |, |, |, 

" gUx,Q 2 )dx = . (5.1) 



o 



Now we can ask what happens to the first moment of the virtual photon structure 
function g1{x,Q 2 ,P 2 ). This can be studied by taking the n — > 1 limit of (|3.16|). 



Note that we have the following eigenvalues of the one-loop anomalous dimension 



matrix 7° =1 : 



X1 =1 = 0, XT 1 = -2/3o, X n N = s =0. (5.2) 



Physically speaking, the zero eigenvalues A+ =1 = A^. 1 = correspond to the 
conservation of the axialvector current at one-loop order, which has the counterpart 
for the unpolarized structure function F 2 ; the conservation of energy momentum 
tensor, A™ =2 = 0. The other eigenvalue A™ =1 = — 2/3q, which is negative, is rather 
an artifact of continuation of the anomalous dimension of the gluon operators to 
n=l, since there is no twist- 2 gluon operator exists for n — 1, in the RG scheme in 
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which only gauge-invariant operators are allowed. But n = 1 gluon operator exists 
in the so called Adler-Bardeen scheme [[3^, [39]. In fact, in the QCD improved parton 
model approach, there is no reason why the n = 1 moment of the polarized gluon 
distribution should not be considered @ . 



In the coefficients 



P 'TTWW ' A -' e " (5.3) 



the special points (|5.2[) would develop the singularities at n = 1, since in those 
coefficients there exist the factors 

1 1 1 (5.4) 



A+' l + A!L/2/V 

Now if we take the limit of n going to 1, we have 

A\ -> finite, A n _ -> 0, A n NS -> finite 

£™ -> 0, B!! -> finite, -> . (5.5) 

However, A\, A^ s , and B™ are multiplied by the following vanishing factors 

(5.6) 

respectively, and thus the terms proportional to P™, Af, and B" in Eq.( p,16|) all 
vanish in the n = 1 limit. Note that these vanishing factors are specific to the case 
of the virtual photon target, and that such factors do not appear when the target 
is real photon. 

Thus the first three terms in the 1st moment vanish irrespective of the RG 
scheme. So we get 

/ dxgj(x, Q\ P 2 ) = ;^C n=1 + 0(a s ) (5.7) 

JO 47T ZfjQ 

Now let us consider C n=1 , which is given by 

C n=1 = 12/3 < e 4 > {B n G + Ala) | n=1 (5.8) 
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As we have seen, the combination (B G + A^ G ) is renormalization-scheme indepen- 
dent P5]| . The results in the MS scheme P, p|, p9 



are 



r>n=l 
B G 



0. 



A 



n=l G 



-2n f 



(5.9) 



The same results have been obtained by Kodaira |28j in the framework of OPE and 
RG method. He set B G =l = 0, observing that there is no gauge-invariant twist-2 
gluon operator for n = 1 and obtained A^ =1 G = —2nf from the Adler-Bell-Jackiw 
anomaly. In the end, we have for the sum rule of the virtual photon structure 
function gj, 



r 1 3« f 

/ ^(x,Q 2 ,P 2 ) = --E4 + 0(« 

JO 71 i=1 



(5.10) 
(5.11) 



where 

Q 2 > P 2 > m^,A 2 , i = l, •••,7i / 

with rrii, the mass of i-ih flavor quark, and rif, the number of active flavors. 

Now it should be pointed out that we can further pursue the QCD corrections 
of order a s to the first moment of g[. In the above equation for the first moment, 
the leading order is 0(1) not of order l/a s (Q 2 ), which is the case for the general 
moments. So we now go to the order a s QCD correction. 



First we take the renormalization scheme of Kodaira [2^]. We write down the 
first moment of gj (x, Q 2 , P 2 ): 

f Q dxgl{x,Q 2 ,P 2 ) = Ci{Q 2 /P 2 rg{P 2 ),a)( 1 { P )\Rt^ 2 = P 2 )\ 1 { P )) 

+ Cf s (Q 2 /P 2 ,^(P 2 ),a)( 7 (p)|i?r(/i 2 = P 2 )\l{p)) (5.12) 

Here it should be emphasized that because of the absence of the gauge-invariant 
n — 1 gluon and photon operators, the mixing problem becomes much simpler. The 
coefficient functions can be given by 



Ct(Q 2 /P 2 ,g(P 2 )-, 

c? s (Q 2 /P 2 ,g(P 2 ) 



a) 
, a) 



= Texp 

where 71(g) is a 2 x 2 diagonal matrix: 

7i(0) : 



(p2) ,yW) / Cf(l,g(Q%a) 

- 9 m 9 M \c? 8 {\,gm,<*) 











iNs{g) 



(5.13) 



(5.14) 
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Here anomalous dimensions are expanded in powers of the coupling constant: 

7(9)=7 <0) I f^ + 7 (1) ( I ^) 2 + O(/), (5.15) 



where 

7$ = 7^ = 0, 7 # = 24C F T, = 24 • *tj± • ^ = 16n,, 7 « = , (5.16) 
and the coefficient functions are 

Cf(l,s(Q 2 ),«) = < e 2 > (1 - =< e 2 > (1 ^ 



Texp f 9(P2) = 1 - —^[?(P 2 ) - g\Q 2 )] • (5.18) 



4 7T 7T 

^(UWV) = 1.(1-^^1) = 1-^1 (5.17) 

4 7T 7T 

r*^ Mg')_ 1 1 7 (1) 

W 2 ) P W) " 16tt 2 2/? ' 
Here we have the finite matrix element of the quark operators between the virtual 
photon states: 

(l(p)\K=i(^ = p 2 )\i{p)) = ^AU, (5-19) 

where 

A n = 6(< e 2 >, 0, < e 4 > - < e 2 > 2 )i^ G . (5.20) 

Now we recall Kodaira's statement that the bare Green's function for the n — 1 
case does not receive divergent corrections but the finite correction connected with 
Adler-Bell-Jackiw anomaly: 

iti G = —4T(R) = -2n f . (5.21) 
Putting all these equations together, we finally obtain the 0(a s ) QCD correction: 



/ dxgj(x, Q 2 , P 
Jo 



2 r>2, 3a 



71 



E4(i as{Q 



i=i \ n / 

2^ 



2 ( g e2 ^ 2 ^s(P 2 ) a s (Q 



+ 0(a 2 s ). (5.22) 



This result is perfectly in agreement with the one obtained by Narison, Shore and 
Veneziano in ref. ||, apart from the overall sign for the definition of gj. 
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Now we show that we obtain the same result for the first moment of gj(x, Q 2 , P 2 ) 
in the MS scheme. Although there exist no gauge-invariant twist-2 gluon- and 
photon-operators for n = 1, the MS calculation of the anomalous dimensions gives 
non-zero results for 'Jq^ 1 and Jq^ 1 - Thus the MS-scheme results rather correspond 
to the QCD parton model approach where the first moments of the gluon and photon 
distributions are as well defined as the other distributions. 

Including the gluon and photon operators, let us start with Eq.(3.13) for the first 
moment of gj(x, Q 2 , P 2 ) in MS scheme: 

j^dxgJ{x,Q 2 ,P 2 ) = ^-A 1 .M 1 (Q 2 /P 2 ,g(P 2 )).C 1 (l,g(Q 2 )) 



+ X 1 (Q 2 /P 2 ,g(P 2 
We expand Ai in powers of g 2 (P 2 ) as 

g\P 2 



C^giQ^ + Cl (5.23) 



Ai = A? + 



16vr 2 



-A 



(i) 



where Ai is given in Eq.(3.23). Using the MS results for the n 
dimensions, we obtain up to the order 0(g 2 ) 

/I- 1 



(5.24) 
1 anomalous 



M 1 



V 



M iSP l9 2 (P 2 ) ~ 9 2 {Q 2 )) x 24C F T f , 
0, 
0. 



\ 


1 / 



(5.25) 



where the second column is irrelevant since Bq 



in MS scheme and thus Cp (1, g) 
starts in 0(a 2 ). Now it is easy to see that the first term of Eq.( [5.23D , to be more 
specific, (^A^MiCi(l, g(Q 2 ))) gives the same result as in Eq. fl5.22p . 

Let us now consider the contributions of other terms. If A^' and A^ s ^ in the 
second term in Eq. (|5.24 ) remain non-zero, then they give the 0(g 2 ) contribution. 
But Ai = A^ s ^ = due to the non-renormalization theorem [[|l|] for the triangle 
anomaly, so its contribution is at most in 0(a 2 ). The contribution of the second 



term in Eq. ( |5.23|) is also in 0(a 2 ), since K\ = K± = and we expect 



K 



(2),n=l 
NS 



o , 



(5.26) 



for the three-loop mixing anomalous dimensions which are implied from the fact 
that the three-loop ^c™ -1 



4C 



20 



Finally we expand the third term of Eg. ( |5.23| ), C*7, as 

C7(l, g, a) = ^ [B^ + ^X ] ' n=l + ■■] ( 5 - 27 ) 

where B^' n=1 = 5" =1 in Eq.(3.15). We already know that B^ n=1 = in MS 
scheme. On the other hand, the two- loop (0(a 2 )) coefficient function for the polar- 
ized gluon has been caluculated in the MS scheme by Zijlstra and van Neerven |30| . 
It is made up of two terms, one proportional to factor CpTfUf and the other pro- 
portional to factor Cj^TfUf. The first moments of both terms turn out to vanish. 
It can be shown that the two-loop (0(aa s )) coefficient function for the polarized 
photon, flW, is obtained from the two- loop gluon coefficient function, by picking 
up the term with the CfTjUj factor and by modifying the group factors. Thus we 
conclude that the first moment of B^ is zero. In the end, C\ does not give 0(1) 
nor 0(a s ) contributions to the first moment of gj (x, Q 2 , P 2 ). This means that we 
arrive at the same result for the 1st moment of gj given in ( |5.22| ) in the MS scheme. 

6 Numerical analysis 



We now perform the inverse Mellin transform of the equation ( |3.16| ) to get g^ as a 
function of x. The n-th moment is denoted as 

M(n,Q 2 ,P 2 ) = [ x n - 1 gJ(x,Q 2 ,P 2 )dx. (6.1) 
Jo 

Then by inverting the moments ( |6.1| ) we get 

1 rC+ioo 

gj(x, Q 2 , P 2 ) = — M(n, Q 2 , P 2 )x~ n dn (6.2) 

In order to have better convergence of the numerical intergration, we change the 
contour in the complex n-plane from the vertical line connecting C—ioo with C -Moo 
(C is an appropriate positive constant), introducing a small positive constant e, to 

n = C — e\y\ + iy, — oo < y < oo (6.3) 

Hence we have 

gUx,Q\P 2 ) = ^- T M{C-ey + iy,Q 2 ,P 2 )e-^y^°^\i~e)dy 
Lix% Jo 
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+ -L [° M(C + ey + iy,Q 2 ,P 2 )e- {c+£y+iy)lo ^ x \t + e)dy 

2lTl J-oo 

1 r°° r 

= - / \Re{M(z,Q 2 ,P 2 )e~ zlosix) } 

71 JO 

— elm{M(z,Q 2 , P 2 )e- zlog{x) }] dy 
z = C — ey + iy . (6.4) 

In Fig. 3 we have plotted, as an illustration, the result for n/ = 3, Q 2 = 30 GeV 2 
and P 2 = 1 GeV 2 for the QCD scale parameter A = 0.2 GeV. The vertical axis 
corresponds to 

gl{x ) Q 2 1 P 2 )/—n f <e i >\n Q2 



TT f P 2 
Here we have shown four cases; the Box (tree) diagram contribution; 



7 (Box) 

9i 1 



x,Q 2 ,P 2 



(2x-l)—n f < e 4 > In — , 



TT 



(6.5) 



(6.6) 



the Box diagram contribution including non-leading correction ignoring quark mass 



9i 



7(Box,non— leading) 



x,Q 2 ,P 2 



3a 



TT 



-n f < e > 



O 2 

(2x - 1) In ^- - 2(2x - l)(\nx + 1) 



(6.7) 

the leading-order (LO) QCD correction and the next-to-leading order (NLO) QCD 
correction. We observe that the NLO QCD correction is significant at large x as 
well as at low x. We have also studied other examples with different Q 2 and P 2 . 
In Fig. 4 we have plotted the case for Q 2 = 100 GeV 2 with P 2 = 1 GeV 2 . Another 
case for Q 2 = 30 GeV 2 with P 2 = 3 GeV 2 is shown in Fig. 5. We have not seen 
any sizable change for the normalized structure function ( |6.5| ) for these different 
values of Q 2 and P 2 . We examined the rif — 4 case as well. It is observed that the 
normalized structure function is insensitive to the number of active flavors. Here we 
have not directly taken into account the heavy quark mass dependence, but rather 
confined ourselves to the above kinematical region. It turns out from the numerical 
analyis as well as from theoretical arguments that, as P 2 increases, the NLO QCD 
result approaches the Box contribution including the non-leading correction, as in 



the unpolarized structure function |TT|, [12 . 
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Now let us cosider the real photon case P 2 = 0. The structure function can be 
decomposed as 

gj(x, Q 2 ) = gl(x, Q 2 ) | pert . + gl(x, Q 2 ) | n0 n-pert. (6.8) 

The second term can only be computed by some non-perturbative method like lat- 
tice QCD, or estimated by vector meson dominance model (VMD). The first term, 
the point-like piece, can be calculated in a perturbative method. Actually, it can 
formally be recovered in our analysis by setting P 2 = A 2 in Eq.( [3.16|) . In Fig. 6, 
we have plotted the point-like piece of g\ of the real photon. The LO QCD result 
coincides with the previous result obtained by Sasaki in The NLO result is 

qualitatively consistent with the analysis by Stratmann and Vogelsang | i~0fl . Finally, 
a comment on the n = 1 limit of the real photon structure function is in order. In 
the case of the unpolarized structure function i 7 ^ we have a singularity of A n _ at 
n = 2 due to the vanishing of A™ at n = 2 which leads to the negative structure 
function | 18fl . As discussed in refs. [[42] , 43 , [44] we have to introduce some regulariza- 



tion prescription to recover positive structure function. For the polarized case we 
do not have such complication at n = 1 as we have seen in sect. 5. 



7 Conclusion 

Here in the present paper, we have investigated virtual photon's spin structure 
function gJ(x,Q 2 ,P 2 ) for the kinematical region A 2 C i 52 < Q 2 , in the next-to- 
leading order in QCD. We presented our arguments both in the framework of OPE 
supplemented by RG method and in the DGLAP equation approach. The results 
are shown to be independent of the renormalization scheme. 

The first moment of gj for the virtual photon is non- vanishing in contrast to the 
vanishing first moment for the real photon case. We can go a step further to the 
order a s which has turned out to reproduce the previous result of Narison, Shore 
and Veneziano ||, and the result is RG scheme- independent. 

The numerial evaluation of gj by the inverse Mellin transform was performed. 
The result shows that the NLO QCD corrections are significant at large x and also 
at small x. The numerical analysis can also be applied to the point-like component 
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of the real photon structure function. The result is qualitatively consistent with the 
previous analysis. 

Although we have neglected in our kinematical region, we should also consider 
the power corrections of the form (P 2 /Q 2 ) k (k = 1, 2, • • •), which are arising from the 
target mass effects as well as from higher-twist effects. 

In the present paper we only presented the result for the polarized photon struc- 
ture function gj itself. In the course of the parton model analysis, we also obtain 
the polarized parton distributions [0, [45[] of the longitudinally polarized photon, for 
the case of virtual photon, which will be discussed elsewhere. 

As a future subject, it would be intriguing to study another spin structure func- 
tion gl which only exists for off-shell photon. In the OPE language, the twist-2 as 
well as twist-3 operators contribute to the QCD effects for gl , which are now under 
investigation. 
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Appendix 



A Notation for anomalous dimensions 

To the lowest order in a, the anomalous dimension matrix has the form 

ln(g,a)=( Ji 9 2 2) ° n ) (A.l) 



K n (g 2 ,a) 0, 

where 7 n (g 2 ) is the usual 3x3 anomalous dimension matrix in the hadronic sector 

(1U9) 73* (0) \ 
%{g) = iUg) laaig) o , (a.2) 
V r N s(g)J 

and K n (g,a) is the three-component row vector 

K n (g,a) = (K^(g,a), K%(g,a), K n NS {g,a)) (A.3) 

representing the mixing between photon and three hadronic operators. 
The anomalous dimensions are expanded as 

Ug) = ^ + jS^^ ) + 0{9% (a - 4) 

K nM = -^ K °--J^y K ^+0(eY). (A.5) 

The one-loop anomalous dimension matrix 7° can be expressed in terms of its eigen- 
values A" (i — +, — , NS) as 

f n = E kpi 1 (a.6) 

i=+ ,-,NS 

where 

\n 1 f 0,n 0,n 1 \/ 0,n 0,n\2 , a 0,ra 0;n\ l / 2 \ / A m 

A ± = 2l 7 ^ + ^gg± -7gg) + 4 7<k?7gVJ } ( A -7) 

A^5 = 7^ (A.8) 
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and P™ are the corresponding projection operators, 

1 / liptp ~~ ^ 



pn _ 

± \% - A™ 



V 



0,n 





0,n 

7g'S-A" 



0\ 







pn 






1 



(A.9) 
(A.10) 



B Explicit expressions for anomalous dimensions 
B.l One- loop order 



where 



and 



0,n 

7gV 



0,n 

Igg 



-877, 
-4C F 



n(n + 1) 



+ 451(71) 



2C7 



n-l 
f n(n + 1) ' 
n + 2 
n(n + 1) ' 
11 £ 



3 n(n + 1) 



+ 451(71) 



+ -Tt 

3 ' 



5i(n) = E^. 



(B.l) 
(B.2) 
(B.3) 
(B.4) 

(B.5) 



(7^ = 3, C F = -, T f = T R n f =n f /2 



with n/ being the number of flavors. 



(B.6) 



J U > NS 



= 24n f < e 2 > 



n — 1 
n(n + 1) 



0,n 
TVS 



24 n/(< e 4 > - < e 2 > 2 ) 



n — 1 
n(n+ 1) 



(B.7) 
(B.8) 

(B.9) 



26 



B.2 Two-loop order g, g| 
B.2.1 Non-singlet sector 

Ins" = 8C| A n NS + 8C A C F B n NS + 8C F T f D n NS , (B. 10) 



with 
.4 



B 



where 



' vs 8 n n + 1 n 2 (n + l) 2 n 3 (n + 1) 3 

r_4 4 _2_ 2 2 2 

+ ^ j l n + n + 1 + n 2 + (n + l) 2 n 3 + (n + l) 



3 , 



+£i(n)(4 - T-T^) +^(n)(3 - - + + 4^(n 
v n^ (n + l)^ 7 v n n + 1 

T7 /2 2 \ ?? 

+^(-)(- - - 4Si(n)) - S' 3 (-) + 8S(n) (B.ll) 

17 187 1 187 1 17 1 5 1 1 1 



NS ~24 ~18n 18" n + 1 If ^ ~ 6 (n + l) 2 (n + l) 3 

f 2 2 1 1 1 1 | 

+ ^ J U n + 1 n 2 (n + l) 2 + n 3 (n + 1) 3 J 

+^S l (n) + 5 2 (n)(-H + ^ _ 2 _ 
9 v 3 n n + 1 

+ S ^){- 1 - + 7^1 + 25i(n)) + \S'^) - AS{n) (B42) 
1 22 1 22 1 2 1 2 1 

£) n — I I 

^ 6 9 n 9 n + 1 3n 2 3 (n + 1) 2 

-ySi(ri) + ls a (n) (B.13) 



and 



S a (n) = E3. 53 H = E^^ ^H = E i -^^0'), (B.14) 

3=1 ^ 3=1 J 3=1 ^ 



n ! + (-!)» w , l-(-l)" w -l, 

^2(2) = 2 S ^ 2 2" ' ^ ' 

n ! + (-!)* n 1 - (-1)" n - 1 

^3(3) = 2 3 ^ 2 2~ ^ ' ^ ^ 
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B.2.2 Singlet sector 

[1] 7VV : 

with 



lw> — Ins + Tps,^ 



n+1 n 2 (n + 1) 2 n 3 (n+l) E 



[2] 7^g: 



with 



U ^G 



rnn 



-n n + 
22 27 9 

+ TT _ - 

n n+1 n z 



24 22 2 

+ 7 + ~T + 



(n+1) 2 
24 



n + 
2 

— + 

4 

+ -7 + 



n + 1 



(n+1) 3 
24 



n ' n + 1 ' n 2 (n + 1) 2 n 3 (n + l) 3 



-SUn 



n n + 1 (n + 1) 2 



n n + 1 



+5 2 (n)(-- 1 



4 



n n + V 2 ; V n n + 1 



+ 



1 + (-i) B ] d-^Ti) {- 2S *v> + ^ ) + C(2)) 



[3] 1g^- 



with 



An 

A Gip 



IgT = % c f A Gip + 8CUC F B g ^ + 8CfT/ 



HI 

In 



n + 1 n 2 2 (n + 1) 2 



n J 



(n + 1) J 



-S 1 (n)( ~ + ^— + + S!(P)(- - 

v 7 V n n + \ (n + 1) 2 / 1V 7V n n + 1' 
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+S 2 (n)(- - 



B G^ 



41 1 



35 1 
In+1 
10 1 



38 



1 



6 



—S 2 (ri 



3 (n + 1) 2 
4 



(n + 1) ; 



+ 



Jn + 3nTI + ^)~ bM ^n~ nTI^ 

^n n + 1/ 2V 2 ; W n + 1' 
l + (-ir](^-^)(25 2 (n)-^)-C(2)) 



[4] 7gg: 



with 



U GG 



^GG 



161 4 1 4 1 

¥n + 9n~+T + 3 (n + 1) 2 ~ V 



Igg 1 = 8C F T f Dq G + 8C A T f E^ G + %C\ Fq G , 



81 4 1 x 
3n 3n + \> 



_ ( 10 10 10 2 

n n+1 n 2 (n + 1) 2 



4 1 



4 

+ — + 
1 



4 76 1 _76^ 

3 + in ~ ~3^~ 3 (n + l) 2 



(n + l) 3 
20 Q 



(B.23) 



(B.24) 
(B.25) 

(B.26) 

(B.27) 
(B.28) 



^GG - 



67 



1 



8 97 1 97 1 29 1 
3 ~ 18 n + 18n + 1 + T n 2 ~ T (n + l) 2 
67 8 8 



24 



rr 



(n + l) f - 



+S 1 (n)( 



9 



n z 



[n + 1) 2 
2 

T 



)-kS+4S(n) 



+ 



1 + (-!)' 



85 2 (n)( 



1 



1 



n + 1 



25 3 (n) - 45!(n)5 2 (n 
2 



^n 

+25^)(5 l( n)-^ + — 
+C(2)(25 1 (n)-1 + -1 T )-C(3) 



^(|)-45(n) 



(B.29) 
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B.2.3 



K = (K' n i K o n > K N n s) , (B.30) 

with 

KY = -3 n f < e 2 > C F 8D^ G , (B.31) 

K]f = -3 n f < e 2 > C F 8(D GG - 1) , (B.32) 

K]f s = -3n/(<e 4 >-< e 2 > 2 )C F 8D^ G . (B.33) 

B.2.4 Anomalous dimensions at n = 1 (MS scheme) 

7^ =1 = 7j; =1 = (B.34) 

1% =1 = 0, (B.35) 

= ~QC F (B.36) 
99 8 

1°gg =1 = - Y C A + -T f = -2P (B.37) 

^™ =1 = Kl' n=1 = (B.38) 

7 S n=1 = (B.39) 



{1) ' n_1 = 24C F T/ (B.40) 

7$T =1 = (B.41) 

7^ n=1 = + (B.42) 

7 ^ n=1 = 8C F T f + fc A T f - fc 2 A = -2A (B.43) 



^(l),n=l = ^g),n=l = ^W,«=l = Q ( R44 ) 
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C Coefficient Functions 

C n (CX) is the coefficient function of the hadronic (photon) operators |17 

C n (l,g(Q 2 ))= ti^BZ 

\5 N s(i + S W B ns)) 



C n (l,g(Q 2 ),a) 



and 



m = (2/n/)U3 



C.l MS scheme 



n ip — n NS 



2 2 
- + — 

n n 



n-l 



— + 3 )S 1 (n-l)+A y £-S l {j)-AS 2 { 
+ 1 / pi J 



6 

+ - -9 

n 



B n G = 4T f 



-(5' 1 (n) + l)-- + 



n(n + 1) 



n 2 n(n + 1) 



C.2 Momentum subtraction |27|, |28fl 



-2-2 



4 



n n + 1 n 2 (n + l) 2 
+35i(n) -8£ 2 (n) 



5£ 



8T, 



G 



1 2 

n n + 1 
= 



(n + l) 1 



n > 3 
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D Tensor decomposition of virtual photon-photon 
amplitude 



After using parity invariance, time-reversal invariance, and gauge invariance, Brown 



and Muzinich |2(J have shown that there are eight independent tensors, in other 
words, eight-invariant amplitudes for virtual photon-photon scattering. Those eight 
independent tensors, which are free from kinematic singularities and kinematic zeros, 
are given in Eqs.(A3)- (A10) of Ref.f2(J. 



Using these tensors {Iij^vpr > the absorptive part of the forward virtual photon- 
photon scattering amplitude W^pr is decomposed as 

8 

Wp Upr = Y.^^uprAiiw, t h t 2 ) (D.l) 
i=l 

where the Ai are the invariant amplitudes and 

w = p-q, t l = q 2 = -Q 2 ) t 2 =p 2 = -P 2 . (D.2) 

In order to implement crossing symmetry under q — > —q and p, «-> v, we form the 
even combinations, I\, I 2 + 1 3 , I4, I5, I7 + 1$, and the odd combinations, I 2 — I3, 
If — Ig, I' 6 = 2/ 6 — 3wl 7 — wl 8 + {tit 2 /w)(I 2 — J 3 ). It is noted that the combinations 
I 2 — J3 and J7 — Is are anti-symmetric under the interchange of p «-> v and p «-> r, 
while the rest of the combinations are symmetric. In terms of these crossing-even 
and -odd combinations, the amplitude W = J2i=i hA+ is rearranged as follows: 

8 1 

Y^hAi = I 1 A 1 + -{I 2 + h){A 2 + A^) + hA A + hA b 



+ ~(I r + I 8 )(A 7 + A s + 2wA 6 ) + ~I^A 6 

+ ~(I 2 - h){A 2 - A 3 - —A 6 ) + ~(/ 7 - J 8 )(A 7 - A 8 + wA 6 )(B.3) 
2 w 2 

Now #7 = 2W4 is written in terms of invariant amplitudes as 

gj oc aim - ai-n-i (D.4) 
= w 2 {A 2 -A 3 - t ^A 6 )-t 1 t 2 (A 7 -A 8 + wA 6 ) (D.5) 
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which is obtained from Eq.(Al4) of Ref . ||20|| . Here ami (ai-11-1) represents the 
s-channel helicity amplitude for (+1)7+ (±1)7 — ► (+1)7+ (±1)7 • It is noted that 
(A2 — A3 — (tit2/w)Ao) is the invariant amplitude associated with (I2 — I3) , while 
(A? — A 8 + wAq) is associated with (I 7 — J 8 ). In the limit t 2 = p 2 = or in the case 
that w = p ■ q,t 1 = q 2 ^> t 2 = p 2 , the second term t 1 t 2 (A 7 — ^8 + wA 6 ) does not 
contribute to gj. 

In fact we observe that the tensor (I2 — I3) = I is associated to g\ while 
(I 7 — I 8 ) = J_ is associated to g^ ■ It can be shown that 

e^A CT gV CT % = — /- , (D.6) 
p-q 

and in the limit of — q 2 ,p ■ q ^> —p 2 

e m \*q v q a - ^ va \<y%(f - e^Aa £ P T al3 P(3P X = J- ■ (D.7) 
Now using an identity 

+ 9^ + g^afrv (D.8) 



we get 



e^Aag A (P • q Cpr^Pp ~ e p ra(3P l3 p a q a ) 
tfiaXaquq ti/aXcrq/iq t/ivXcrq 



(D.9) 



Hence we have from Eq.(2.3) 



W, upT = _— [(J.)^ gj - {J_)^ pT gl]. (D.10) 

Finally it is interesting to see the relation between the polarized photon struc- 
ture functions gj and g\ and polarized nucleon structure functions g\ and g^. By 
introducing the polarization vectors, e* p and e T for the target photon just like those 



for the gluon target discussed by Gabrieli and Ridolfi we have 

1 

■~2 [ 



W^i =e pT ^p lSs ) . (D.li; 

2 vb 2 l 
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where s is the longitudinal spin vector for the target photon. After using the relation 
p ■ s = 0, we get 



W, 



[IV 



p-q 



e^\ a q X s a g1 + e livXa q x {p ■ qs a - q ■ sp a ] 



7 
92 

p-q 



(D.12) 



which, apart from the factor y \p 2 \, has exactly the same form as Eq.(2.4) of Kodaira 
et al. |27j which defines the polarized nucleon structure functions g± and g?,-, and also 



as Eq.(9) in Ref. [@. 



34 



References 

[1] D. Barber, in proceedings of the "Zeuthen Workshop on the Prospects of Spin 
Physics at HERA", DESY 95-200, p.76, eds. J. B16mlein and W. D. Nowak. 

[2] M. Stratmann and W. Vogelsang, Z. Phys. C74 (1997) 641. 

[3] S. D. Bass, Int. J. Mod. Phys. A 7 (1992) 6039. 

[4] A. V. Efremov and O. V. Teryaev, Phys. Lett. B240 (1990) 200. 

[5] S. Narison, G. M. Shore and G. Veneziano, Nucl. Phys. B391 (1993) 69; 
G. M. Shore and G. Veneziano, Mod. Phys. Lett. A8 (1993) 373; 
G. M. Shore and G. Veneziano, Nucl. Phys. B381 (1992) 23. 

[6] A. Freund and L. M. Sehgal, Phys. Lett. B341 (1994) 90. 

[7] S. D. Bass, S. J. Brodsky and I. Schmidt, Phys. Lett. B437 (1998) 424. 

[8] R. Mertig and W. L. van Neerven, Z. Phys. C70 (1996) 637. 

[9] W. Vogelsang, Phys. Rev. D54 (1996) 2023; Nucl. Phys. B475 (1996) 47. 
[10] M. Stratmann and W. Vogelsang, Phys. Lett. B386 (1996) 370. 
[11] T. Uematsu and T. F. Walsh, Phys. Lett.lOlB (1981) 263. 
[12] T. Uematsu and T. F. Walsh, Nucl. Phys. B199 (1982) 93. 
[13] G. Rossi, Phys. Rev. D29 (1984) 852. 

[14] M. Drees and R. M. Godbole, Phys. Rev. D50 (1994) 3124. 

[15] M. Gliick, E. Reya and M. Stratmann, Phys. Rev. D51 (1995) 3220; Phys. Rev. 
D54 (1996) 5515. 

[16] P. Mathews and V. Ravindran, Int. J. Mod. Phys All(1996) 2783. 
[17] W. A. Bardeen and A. J. Buras, Phys. Rev. D20 (1979) 166. 

35 



[18] D. W. Duke and J. F. Owens, Phys. Rev D22 (1980) 2280. 

[19] E. Witten, Nucl. Phys. B120 (1977) 189. 

[20] R. W. Brown and I. J. Muzinich, Phys. Rev. D4 (1971) 1496. 

[21] M. A. Ahmed and G. G. Ross, Phys. Lett. B59 (1975) 369. 

[22] K. Sasaki, Phys. Rev. D22 (1980) 2143; Prog. Theor. Phys. Suppl. 77 (1983) 
197. 

[23] A. V. Manohar, Phys. Lett. B219 (1989) 357. 

[24] K. Sasaki and T. Uematsu, "Polarized Virtual Photon Structure Function in 
the Next-to-leading Order in QCD ", preprint, [hep-ph/9812380| . 

[25] W. A. Bardeen, A. J. Buras, D. W. Duke and T. Muta, Phys. Rev. D18 (1978) 
3998. 

[26] J. Kodaira, S. Matsuda, T. Muta, K. Sasaki and T. Uematsu, Phys. Rev. D20 
(1979) 627. 

[27] J. Kodaira, S. Matsuda, K. Sasaki and T. Uematsu, Nucl. Phys. B159 (1979) 
99. 

[28] J. Kodaira, Nucl. Phys. B165 (1980) 129. 

[29] G. T. Bodwin and J. Qiu, Phys. Rev D41 (1990) 2755. 

[30] E. B. Zijlstra and W. L. van Neerven, Nucl. Phys. B417 (1994) 61. 

[31] W. A. Bardeen and A. J. Buras, Phys. Lett. 86B (1979) 61. 

[32] A. J. Buras, Rev. Mod. Phys. 52 (1980) 199. 

[33] M. Gliick and E. Reya, Phys. Rev. D25 (1982) 1211. 

[34] G. Altarelli, Phys. Rep.81 (1982) 1. 



36 



[35] W. Furmanski and R. Petronzio, Z. Phys. Cll (1982) 293. 

[36] M. Gliick and E. Reya, Phys. Rev. D28 (1983) 2749; 

M. Gliick, E. Reya and A. Vogt, Phys. Rev. D45 (1992) 3986. 

[37] M. Fontannaz and E. Pilon, Phys. Rev. D45 (1992) 382. 

[38] G. Altarelli and G. G. Ross, Phys. Lett. B212 (1988) 39. 

[39] R. D. Ball, S. Forte and G. Ridolfi, Phys. Lett. B378 (1996) 255. 

[40] G. Altarelli and B. Lampe, Z. Phys. C47 (1990) 315. 

[41] S. L. Adler and W. A. Bardeen, Phys. Rev. 182 (1969) 1517. 

[42] W. A. Bardeen, Proc. 10-th Int. Symposium on Lepton and Photon Interaction 
at High Energies, (Bonn, 1981), ed. W. Pfeil (Bonn University, 1981) p. 432. 

[43] I. Antoniadis and G. Grunberg, Nucl. Phys. B213 (1983) 445. 

[44] D. W. Duke, Proc. XIV Int. Symposium on Multiparticle Dynamics, (Lake 
Tahoe, 1983) eds. P. Yager and J. F. Gunion (World Scientific,1984) p.124. 

[45] M. Gliick, M. Stratmann and W. Vogelsang, Phys. Lett. B337 (1994) 373. 

[46] A. Gabrieli and G. Ridolfi, Phys. Lett. B417 (1998) 369. 



37 



Figure Captions 

Fig. 1 Deep inelastic scattering on a polarized virtual photon in polarized e + e~ 
collision, e + e~ — > e + e~+ hadrons (quarks and gluons). The arrows indicate 
the polarizations of the e + , e~ and virtual photons. The mass squared of the 
"probe "("target ") photon is -Q 2 (-P 2 ) (A 2 < P 2 < Q 2 ). 

Fig. 2 Forward scattering of a virtual photon with momentum q and another virtual 
photon with momentum p. The Lorentz indices are denoted by /i, u, p, r. 

Fig. 3 Polarized virtual photon structure function gJ(x,Q 2 ,P 2 ) to the next-to- 
leading order (NLO) in units of (3cm/ < e 4 > fir) \n{Q 2 /P 2 ) for Q 2 = 30 
GeV 2 , and P 2 = 1 GeV 2 and the QCD scale parameter A = 0.2 GeV with 
rif — 3 (solid line). We also plot the leading order (LO) result (long-dashed 
line), the Box (tree) diagram (2dash-dotted line) and the Box including non- 
leading contribution, Box (NL) (short-dashed line). 

Fig. 4 Virtual photon structure function gJ(x,Q 2 , P 2 ) for Q 2 = 100 GeV 2 , and 
P 2 = 1 GeV 2 with A = 0.2 GeV, n f = 3. 

Fig. 5 Virtual photon structure function gi(x,Q 2 ,P 2 ) for Q 2 = 30 GeV 2 , and 
P 2 = 3 GeV 2 with A = 0.2 GeV, n f = 3. 

Fig. 6 Point-like piece of the real photon structure function gi(x,Q 2 ) in NLO for 
Q 2 = 30 GeV 2 with A = 0.2 GeV, n f = 3 (solid line). Also plotted are the 
LO result (long-dashed line) and the Box (tree) diagram contribution (short- 
dashed line). 
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